Many-body effects on adiabatic passage through Feshbach resonances 
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We theoretically study the dynamics of an adiabatic sweep through a Feshbach resonance, thereby 
converting a degenerate quantum gas of fermionic atoms into a degenerate quantum gas of bosonic 
dimers. Our analysis relies on a zero-temperature mean-field theory which accurately accounts 
for initial molecular quantum fluctuations, triggering the association process. The structure of 
the resulting semiclassical phase-space is investigated, highlighting the dynamical instability of the 
system towards association, for sufficiently small detuning from resonance. It is shown that this 
instability significantly modifies the finite-rate efficiency of the sweep, transforming the single-pair 
exponential LZ behavior of the remnent fraction of atoms T on sweep rate a, into a a power law 
dependence as the number of atoms increases. The obtained nonadiabaticity is determined from the 
interplay of characteristic timescales for the motion of adiabatic eigenstates and for fast periodic 
motion around them. Critical slowing-down of these precessions near the instability, leads to the 
power- law dependence. A Linear power-law V oc a, is obtained when the initial molecular fraction is 
smaller than the 1/N quantum fluctuations, and a cubic-root power-law F oc a 1 ^ 3 is attained when 
it is larger. Our mean-field analysis is confirmed by exact calculations, using Fock-space expansions. 
Finally, we fit experimental low temperature Feshbach sweep data with a power-law dependence. 
While the agreement with the experimental data is well within experimental error bars, similar 
accuracy can be obtained with an exponential fit, making additional data highly desirable. 

PACS numbers: 05.30.Fk, 05.30.Jp, 3.75.Kk 



I. INTRODUCTION 



Many of the most exciting experimental achievements 
in ultra-cold atomic physics in recent years have used 
Feshbach resonances 0, 0, S S ■ Not only arc they 
a tool for altering the strength and sign of the interac- 
tion energy of atoms 0, , they also provide a means 
for converting atom pairs into molecules, and vice versa 
!1I10&@,I3[IJ[3EQD1-A magnetic Feshbach 
resonance involves the collisional coupling of free atom 
pairs (the asymptotic limit at large internuclear separa- 
tion, of the incident open channel molecular potential) 
in the presence of a magnetic field, to a bound diatomic 
molecule state (the closed channel) on another electronic 
molecular potential surface. The difference in the mag- 
netic moments of the atoms correlating asymptotically 
at large internuclear distance to the two potential energy 
surfaces, allows the Feshbach resonance to be tuned by 
changing the magnetic field strength 0. Sweeping the 
magnetic field as a function of time, so that the bound 
state on the closed channel passes through threshold for 
the incident open channel from above, can produce bound 
molecules. This technique has proved to be extremely ef- 
fective in converting degenerate ato mic gase s of fermions 
H0,I1 IS EIH El and bosons OElllll into bosonic 
dimcr molecules. Fermions are better candidates for Fes- 
hbach sweep experiments due to the relatively long life- 
times of the resulting bosonic molecules, originating from 
Pauli blocking of atom-molecule and molecule-molecule 
collisions when the constituent atoms are fermions |l6j . 

Here we consider the molecular production efficiency of 
adiabatic Feshbach sweep experiments in Fermi degener- 



ate gases. We determinine the functional dependence of 
the remnent atomic fraction T on the Feshbach sweep rate 
a, following the treatment in a previous Letter ex- 
tending the calculations, and presenting a more detailed 
account of the theoretical methodology. 

The Fermi energy is the smallest energy scale in the 
system in the fermionic Feshbach sweep experiments of 

Refs. II Eli, 03 III [Li. Hence, we treat the fermions 

theoretically as occupying the lowest possible many-body 
state consistent with symmetry considerations arising 
from the method of preparation. Consequently we as- 
sume that the quantum states are filled up to the Fermi 
energy in a fashion consistent with the symmetry prop- 
erties of the gas. In this sense, the gas can be thought of 
as a zero temperature gas. 

The Landau-Zcncr (LZ) model [3 is the paradigm 
for explaining how transitions occur in the collision of 
a single pair of atoms in a Feshbach sweep experiment. 
Theoretical analysis of molecular production efficiency 
in Feshbach sweep experiments in a gas phase have been 
based on Landau-Zener theory £|. Exponential fits have 
also been carried out for experimental molecular effi- 
ciency data. Fig. ^ shows data from experiments on a 
quantm degenerate gas of 6 Li atoms pj, plotting the 
remaining fraction of atoms (red squares) as a func- 
tion of the inverse magnetic sweep rate. The inset of 
Fig. ^ includes an exponential fit (blue dashed curve) , 
T = 0.479 exp(-a/1.3) +0.521, taken from Ref. 0. 
While the exponential curve lies well within all experi- 
mental error bars, the data fits a linear power-law de- 
pendence (green curve) to the same level of accuracy. In 
what follows we provide the theoretical detail required to 
obtain the linear power-law fit in Fig. ^ We show that 
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FIG. 1: (Color online) Fraction of remnant atoms, F, versus 
inverse ramp speed 1/B across the 543 G resonance of 6 Li. 
The experimental data (red squares) of 0, which saturates at 
a remnant fraction of 1/2 [Til , and the mean-field calculations 
(solid green curve) obey a linear dependence on sweep rate 

beyond 0.5 ms/G. is multiplied by 0.5 ms/G to scale the 
abscissa for the calculated results. Also shown as a dashed 
blue line in the inset is the best exponential fit to the data, 
T = 0.479 exp(-a/1.3) + 0.521. 



due to the nonlinearity of the reduced single-particle (i.e. 
mean-field) description of the many-atom system, insta- 
bilities are made possible. These instabilities result in 
the failure of the standard LZ theory when the number 
of atoms is large. We also predict two different power-law 
behaviors, r oc a and T tx a 1 / 3 , depending on the initial 
state of the system prior to the Feshbach sweep. 

The paper is organized as follows. In Section |H] we in- 
troduce the Feshbach model-Hamiltonian and the main 
approximations used. In Sec. IIIII we describe the classi- 
cal phase space obtained for the above model. Section 
^describes how the molecular production efficiency de- 
pends on the stability of the fixed points for the equations 
of motion. In Sec. [V] we describe the role of quantum 
fluctuations which lead to the linear dependence of the 
molecular production efficiency on the sweep rate. The 
analysis of sections llVI and Ivl is verified by the exact nu- 
merical calculations presented in section IVT1 Section lYllI 
contains summary and conclusions. 



II. THE ZERO TEMPERATURE MODEL 
SYSTEM 

Experiments on molecule production in slowly swept, 
broad Feshbach resonance systems H are well ex- 
plained by employing a thermodynamic equilibrium 
model hj. The narrow Li resonance, traversed much 
more rapidly 0, is not expected to fit such a descrip- 
tion. We consider the zero temperature limit to model 
such experiments. At low temperatures one can use a sin- 
gle bosonic mode Hamiltonian [HE1 El El El El El 
because of the Cooper instability which singles out the 
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FIG. 2: (Color online) Many-body collective dynamics of adi- 
abatic passage from a fermionic atomic gas into a molecular 
Bose Einstein Condensate (BEC) for five pairs of fermionic 
atoms, (a) Sweep rate a = 2g 2 N, (b) Sweep rate a = g 2 N/4. 
Overall efficiency is independent of atomic dispersion in both 
(a) and (b). 



zero momentum mode of the molecules produced. Thus, 
we take the Hamiltonian to be 



H 



e * c l,* c *.> 

k,CT 



(1) 



where = h 2 k 2 /2m is the kinetic energy of an atom with 
mass m, and g is the atom-molecule coupling strength. 
The molecular energy £ (t) = Ep — at is linearly swept at 
a rate a > where Ep denotes the fermi energy of the 
atoms, through resonance to induce adiabatic conversion 
of fermi atoms to bose molecules. The annihilation opera- 
tors for the atoms, Ck,<r, obey fermionic anti-commutation 
relations, whereas the molecule annihilation operator bo 
obeys a bosonic commutation relation. 

The Hamiltonian can be further simplified by neglect- 
ing fermionic dispersion. This approximation has been 
commonly used 0, and accounts for the use of a sim- 
ple two-level LZ model, as opposed to a multilevel one, for 
such systems. To justify this assumption we have con- 
ducted many-body numerical simulations to determine 
the effect of fermionic dispersion on the adiabatic conver- 
sion efficiency. Fig. [21 shows exact numerical results for 
the adiabatic conversion of five atom pairs into molecules, 
for different values of the atomic level spacing (and hence 
of the Fermi energy Ep). It demonstrates that the final 
adiabatic conversion efficiency is completely insensitive 
to the details of the atomic dispersion. It is evident that, 
while the exact dynamics depends on Ep, and the levels 
are sequentially crossed as a function of time as the bound 
state crosses the level energies, the same final efficiency 
is reached regardless of the atomic motional timescale 
(i.e., regardless of level spacing). In particular, the figure 
shows that in the limit as a — ► 0, it is possible to convert 
all atom pairs into molecules. This is a unique feature 
of the nonlinear parametric coupling between atoms and 
molecules which should be contrasted with a marginal 



3 



conversion efficiency expected for linear coupling in the 
multi-level LZ model. 

E mployi ng the degenerate model with e k = e for all k 
|24l 12a l2q . it is convenient to rewrite the Hamiltonian in 
terms of the following lowering, raising and z-component 
operators |2fill27| : 



J- 



k c k) TC-k,_t 



(N/2) 3 / 2 



J+ 



k,i C k,T 



b 



(N/2) 3 / 2 



Jz = 



N 



(2) 



where N = 2&q6 + J2k a c ic <r c k,o- is the conserved to- 
tal number of particles. It is important to note that 
J-,J+,J Z do not span SU(2), since the commutator 
[J+,J-] yields a quadratic polynomial in J z (despite 
the fact that the commutators [J+,J Z ] and [J-,J Z ] 
have the right commutation relations). The operators 
Jx = J+ + J- and J y = —i(J+ — J-) can also be de- 
fined. Up to a c-number term, Hamiltonian takes the 
form 



(3) 



where A(t) = 2e — £(t) = at. Defining a rcscaled time 
r = \/N~gt, and assuming a filled Fermi sea (i.e., that 
the number of avialable fermionic states is equal to the 
number of particles), we obtain the Heisenberg equations 
of motion, 



(It 



J x = S(t)J 



y > 



dr v 



dr Jz 



3V2 



-5(t)J x + ^-(J z - 1) 



N V 

V2Jy, 



J Z ) 



(4) 



which depend only on 



A(t)/VNg = (a/g 2 N)i 



the scaled detuning S(t) = 
It is interesting to note that 



exactly these equations of motion are obtained for the 
two-mode atom-molecule BEC [27J where, for the bosonic 
case, lowering, raising and z-component operators are de- 
fined as 



J- 



(N/2) 3 / 2 



J+ 



a\a\bo 
(N/2) 3 ' 2 



Jz = 



N 



(5) 



where a\ and are bosonic annihilation operators obey- 
ing bosonic commutation relations. In these definitions, 
the sign of the operator J z has been reversed relative to 
Eq. (|2J), mapping fermionic association to bosonic disso- 
ciation [2II25II26II. 
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FIG. 3: (Color online) Two dimensional surfaces depicting 
classical phase space: (a) Bloch s phe re of a two-mode Bose- 
Josephson system, as in Refs. [271 . l29l (b) zero single-particle 
entropy surface of the atom-molecule system, as in Eqs. ©. 



III. CLASSICAL PHASE SPACE 

The mean-field limit of Eqs. Q is given by replacing 
Jx,Jy, and J z by their expectation values it, v, and w 
which correspond to the real and imaginary parts of the 
atom-molecule coherence and the atom-molecule popula- 
tion imbalance, respectively. Since quantum fluctuations 
in J z are of order 1/7V, it is also consistent to omit the 
quantum noise term v2(l + Jz)/N in the equation for 
djy/dr in as long as J z is of order unity. For small 
w however, when the molecular population is of the order 
of its quantum fluctuations, this quantum term becomes 
dominant and will have a significant effect on sweep effi- 
ciency, as will be shown in Sec. [3 

In the classical field limit, the equations of motion 



dr 
d 

d 



S(t)v , 

3 ^2 , / 1 

—0(t)u H — (w — 1) I w + - 



—w = V2v , 



dr 



(6) 



depict the motion of a generalized Bloch vector on a two- 
dimensional 'tear-drop' shaped surface, (Fig. Ob), deter- 
mined by the constraint, 



v 2 = ^(w-l) 2 (w + l), 



(7) 



corresponding to the conservation of single-pair atom- 
molecule coherence. The peculiar shape of this equal- 
single-pair-entropy surface is a result of the commutation 
relations for the operators J- , J + , J z . For comparison, 
the two-level spin Hamiltonian may be written only in 
terms of SU(2) generators j23, l29j| and the mean-field 
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FIG. 4: (Color online) The adiabatic eigenvalues of a linear 
LZ problem (a) are contrasted with the adiabatic level scheme 
of the atom-molecule nonlinear system (b). For a critical de- 
tuning, one of the adiabatic eigenvalues in (b) splits into two, 
resulting in the emergence of an additional hyperbolic fixed 
point. 
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motion is restricted to the surface of a Bloch sphere (see 
Fig. EJi) as the mean-field approximation allows for the 
factorization of the ST/ (2) Casimir operator J 2 = J Z (J Z — 
1) + J+ J_ = j(j + 1) into the constraint 



2 i 2 i 2 -2 

u + v + w = j 



(8) 



The surface defined by the constraint Q should be 
viewed as the atom-molecule equivalent of this Bloch 
sphere. Accordingly, we proceed by following the meth- 
ods of Ref. which correspond to a Bloch sphere 
like phase space, generalizing them to the atom-molecule 
parametric coupling case. 

Since the constraint of Eq. (jJJ) restricts the dynamics 
to the two dimensional surface depicted in Fig. 03, it is 
readily seen that in the mean field limit, Hamiltonian @ 
is replaced by the classical form 

H{w,9;S) = (Sw + v/(l + w)(l - w 2 ) cosOj . 

(9) 

Here the Hamiltonian is expressed only in terms of 
the relative phase between atoms and molecules 9 = 
arctan(«/u), corresponding to the azimuthal angle in 
Fig. and the atom-molecule population difference w, 
corresponding to its cylindrical axial coordinate. 

The eigenvalues of the atom-molecule system at any 
given value of 6 correspond to the extrema (wq,9q) of 
the classical Hamiltonian ©, or equivalently, to the fixed 
points (uq, vq, wo) of the mean-field equations JBJl]: 



VQ = 



V2 



(w - 1) (3w + 1) = Su Q ■ (10) 



The number of fixed points depends on the parameter 
5. The point u = v = 0, w a — 1 is stationary for any 
value of 6. Using Eqs. Q and lfTfj]l. other fixed points 
are found from the solutions of 



(3w + 1) : 
4(w + 1) 



(11) 



FIG. 5: (Color online) The ten lowest energy eigenvalues for 
Hamiltonian Q, drawn for the degenerate case, tk = 0. 



in the domain wo € [—1,1]. Consequently, it is evident 
that for \S\> V2 there are only two stationary solutions, 
one of which is wo = 1. However, this stationary point 
bifurcates at the critical detuning of 5 C = \/2, so that 
for \S\ < 5 C , there are three eigenstates, as depicted in 
Fig. |3J). In contrast, in the linear LZ problem (Fig. 
eigenvalue crossings are avoided, and there are only two 
eigenstates throughout. 

The relation between the reduced sigle-particle picture 
of Fig.^Ja and the full many-body sytem it approximately 
represents, is illustrated in Fig. [SJ depicting numerically 
calculated eigenvalues for ten atom pairs as a function of 
(5, when Ep = 0. One can envisage how, when adding 
more and more energy levels, finally collapsing the levels 
to a single curve, the curve structure of Fig. ^emerges. 
The bifurcation of the all-atoms mode is shown to emerge 
from its quasi-degeneracy with slightly higher many-body 
states with a few more molecules. 

Stability analysis of the various fixed points, preformed 
by linearization of the dynamical equations (J3J about 
(ito, vq, wo), yields the frequency Qo of small periodic or- 
bits around them: 



n 



gVN 



V^ 2 + (1 - 3w ) 



(12) 



From Eq. (|12|l it is clear that the characteristic oscilla- 
tion frequency about the stationary point (0, 0, 1) will be 
\/5 2 — 2. Thus, for S > S c the point uq = Vo = 0, Wo = 1 
is an elliptical fixed point, whereas for 5 < 5 C it becomes 
a hyperbolic (unstable) stationary point, with an immag- 
inary perturbation frequency. For the remaining eigen- 
values, we can use (|11H to obtain 



gVN 



1 (1 - wq)(3wq + 5) 
4(w + 1) 



(13) 



5 



giving real f^o for all wq G [—1,1], approaching zero as 
w a — > 1. Consequently, for \S\ < S c there are a total of 
two elliptical fixed points, whereas for |<5| < S c there are 
two elliptical and one hyperbolic fixed points. 

The fixed point analysis carried out so far is summa- 
rized in Fig. where we plot the phase-space trajecto- 
ries, corresponding to equal-energy contours of Hamilto- 
nian JjJJl, for different values of 5. As expected from |jSJ, 
the plots have the symmetry (w, 8; S) <-> (w, 9 + tt; —6). 
For sufficiently large detuning. \5\ > \/2, Eq. Ijlip has 
only one solution in the range — 1 < Wo < 1. There- 
fore there are two (elliptic) fixed points, denoted by a 
red circle corresponding to the solution of Eq. , and 
a blue square at (0,0,1). As the detuning is changed, 
one of these fixed points (red circle) smoothly moves 
from all-molecules towards the atomic mode. At the 
critical detuning S = —\/2 a homoclinic orbit appears 
through the point (0,0,1) which bifurcates into an un- 
stable (hyperbolic) fixed point (black star) remaining on 
the atomic mode, and an elliptic fixed point (blue square) 
which starts moving towards the molecular mode. Con- 
sequently, in the regime \S\ < \/2 there are two elliptic 
fixed points and one hyperbolic fixed point, correspond- 
ing to the unstable all-atoms mode. Another crossing 
occurs at S = \[2 when the fixed point which started 
near the molecular mode (red circle) coalesces with the 
all-atoms mode (black star). Plotting the energies of the 
fixed points as a function of detuning, one obtains the 
adiabatic level scheme of Fig. 0Jd. 

As previously noted, for \S\ < v2 the period of the ho- 
moclinic trajectory beginning at (0,0,1) diverges. This 
divergence significantly affects the efficiency of an adi- 
abatic sweep through resonance because the linear re- 
sponse time to a perturbation in the Hamiltonian be- 
comes infinitely long. Consequently, the sweep is never 
truly adiabatic, and its expected efficiency is lower than 
the corresponding LZ efficiency. This effect is discussed 
in the following section. 



IV. EFFECT OF FIXED-POINT INSTABILITY 
ON ADIABATIC SWEEP EFFICIENCY 



Having characterized the classical phase-space struc- 
ture for the parametrically coupled atom-molecule sys- 
tem, we turn to the process of adiabatically sweeping 
the detuning 6 through resonance, converting fermion 
atoms into bose molecules. As usual, adiabatic follow- 
ing involves two typical timescales: the sweep time scale 
given by the inverse sweep rate 1/a and the timescale 
associated with the period 1/flo of small periodic orbits 
around the fixed point, given in Eqs. (|12f) and l|13f) . Slow 
changes to the Hamiltonian (e.g., by variation of the de- 
tuning S) shift the adiabatic eigenvalues, as depicted in 
Fig. Starting with such an eigenvalue (e.g., the all- 
atoms mode for an initial large negative <5), the state of 
the system will only be able to adiabatically follow the 
fixed point if its precession frequency about it, f^o, is large 



8=-10 




FIG. 6: (Color online) Equal-energy contours of Hamiltonian 
© plotted as a function of w and 8 for different detunings 
S. w = 1 is all atoms and w = —1 is all molecules. The 
various fixed points corresponding to adiabatic eigenvectors 
are marked by (blue) squares, (red) circles and (black) stars. 



compared to the rate at which it moves. The adiabatic 
conversion efficiency is related to the magnitude of this 
precession, which in turn is proportonal to the classical 
action accumulated during the sweep. 

The relation between the sweep conversion efficiency 
and the accumulated classical action is illustrated in 
Fig. [3J Consider first the SU (2) case of Fig. where 
mean- field motion is restricted to the surface of the Bloch 
sphere u 2 + v 2 + w 2 = 1. This illustration applies to the 
standard LZ case Q§] as well as to the nonlinear Bose- 
Josephson system |27l l29l l30j . Having started from the 
'south pole' (0, 0,-1) and carried out the sweep through 
resonance, the classical state Bloch vector (u, v, w) car- 
ries out small precessions about the final adiabatic eigen- 
vector which (for sufficiently large final detuning) is par- 
allel to the w axis. The surface-area enclosed within this 
periodic trajectory is just the action, AI cx u 2 + v 2 , ac- 
cumulated during the sweep. In the extreme limit of 
perfect adiabatic following, the precession approaches a 
point trajectory, having zero action. Larger nonadia- 
baticity leads to larger precession amplitude, and hence 
to larger accumulated action. The remanent fraction in 
the initial state, r = (1 — w), is directly related to AI 
by the conservation rule u 2 + v 2 = (1 + w)(l — w), which 
near w = 1 can be linearized to give AI w 2r. This is 
the expression usually used in calculating LZ transition 
probabilities [Tll29l|. 

For the atom-molecule parametric system of Fig. ^jp, 
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the situation is slightly different, since u 2 + v 2 = (1 + 
w)(l — w) 2 , leading to a square root dependence of the 
remnent atomic fraction T on A J near w = 1. In order 
to estimate A/, we transform w, 8 into action-angle vari- 
ables /, 4>. In terms of these variables the non-adiabatic 
probability T at any finite sweep rate a is given by 



AI 
~2 



R(I, 



(14) 



where R(I, <ft) is related to the generating function of the 
canonical transformation (w,9) — ► {I,4>) psl. l29[. 

Equation l|14|l reflects our discussion on characteristic 
timescales. In order to attain adiabaticity, the rate of 
change of the adiabatic fixed points through the vari- 
ation of the adiabatic parameter A, R(I, 4>) A, should 
be slow with respect to the characteristic precession fre- 
quency <f> = VIq about these stationary vectors. The ac- 
tion increament is proportional to the ratio of these two 
timescales. 

As long as 4> does not vanish, the accumulated action 
can be minimized by decreasing A. For a perfectly adia- 
batic process where A/0 — > 0, the action is an adiabatic 
invariant, so that a zero-action elliptic fixed point evolves 
into a similar point trajectory. For finite sweep rate, the 
LZ prescription [T^j evaluates the integral in l|14|) by in- 
tegration in the complex plain, over the contour of Fig. [7| 
noting that the main contributions will come from sin- 
gular points, where 4> approaches zero and the integrand 
diverges. Since for a linear LZ system there are no insta- 
bilities, all such singularities are guarenteed to lie off the 
real axis, leading to the exponentially small LZ transition 
probabilities. 

The situation changes for nonlinear systems, where in- 
stabilities arrise. In section lHfl we have shown that for the 
atom-molecule system with fermion atoms, the all-atoms 
mode becomes unstable to association when the detun- 
ing hits the critical value of 8 C = \/2. From Eq. Ijl2(l it is 
clear that the characteristic frequency 4> = f2o vanishes 
near u>o = 0. Consequently, there are singular points 
of the integrand in (|14f) lying on the real axis. In what 
follows, we show that these poles on the real axis lead 
to power-law dependence of the transfer efficiency on the 
sweep rate. 

In order to evaluate the integral i|14f> . we need to in- 
vestigate how the characteristic frequency cb = fio de- 
pends on the action-angle <f> near the instability point 
(u = 0, ?; = 0, w = 1), where most action (and hence 
most nonadiabatic correction) is accumulated. It is ev- 
ident from Eq. (|12|l that the precession frequency near 
that point vanishes as 



g^N(l - w ) . (15) 

Differentiating Eq. (|11J) with respect to t, we find that 
the rate of change of adiabatic fixed-points due to a linear 
sweep is, 



w 



4a (w + l) 3/2 
gy/N 3w + 5 



(16) 
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FIG. 7: (Color online) Contour of integration in LZ theory, 
for calculating the integral in Eq. 1141 . All singularities lie off 
the real axis. 



Having found wq, we can now find the explicit form of 
the transformation from wq to the action-angle variable 
(j), near the instability. The action-angle may be written 
as 



cj)dt 



dwo 



Wq 



(17) 



In the vicinity of the singularity (wq = 1) we have S7o 
gy/N(l — wq) and wq ~ V^a/gy/N, resulting in 



g 2 NV2 



(l-^o) 3/2 



(18) 



Thus, as the adiabatic eigenstate approaches wq — 1 the 
angle </> vanishes as (1 — u>o) 3 ^ 2 whereas the characteris- 
tic frequency 4> approaches zero as (1 - wo) 1 /' 2 . Conse- 
quently, we finally find from (|15f) and i|18|) that near the 
singularity, (f> is given in terms of < 



as 




(19) 



Substituting fH^I and A = a into Eq. iffl]) we find that 
the nonadiabatic correction depends on a as 



r oc a 1 / 3 



V. QUANTUM FLUCTUATIONS 



(20) 



So far, we have neglected the effect of quantum fluc- 
tuations, which may be partially accounted for by the 
c-number limit of the source term (\/2/N)(l + J z ) in 
Eqs. As a result, we found in the previous section 
that wq does not vanish as wq approaches unity, and the 
remaining atomic population scales as the cubic root of 
the sweep rate if the initial average molecular fraction is 
larger than the quantum noise. However, starting purely 
with fermion atoms (or with molecules made of bosonic 
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atoms), corresponding to an unstable fixed point of the 
classical phase space, fluctuations will serve to trigger the 
association process and will thus initially dominate the 
conversion dynamics. 

In order to verify that such quantum fluctuations can 
be accurately reproduced by a 'classical' noise term near 
the unstable fixed point w = I, we compare the onset 
of instability from exact many-body calculations to the 
onset of mean-field instability according to the revised 
mean- field equations, 



dr 
d 



dr 



5{t)v, 

—S(t)u 
, V2 



3V2 



(w- 1) \ w+ - 



N 
V2v 



(1 



(21) 



have retained the 0(1/N) noise term 
to). The results, shown in Fig. [HI show ex- 



where we 
(V2/N)(l 

cellent agreement in the early-time dynamics, indicating 
that the mean-field noise term gives the correct behavior 
near the instability point. 

Having established the accuracy of the noise term in 
Eqs. I|21|) we proceed to investigate its effect on sweep 
efficiencies. When this additional term is accounted for, 
Eq. must be replaced by 



6 = 



y/lVQ + 1 



3wq + 1 w + 1 



(22) 



4 iVf> -l), 

This expression reduces to our previous result in Eq. 
when the second term on the r.h.s. of Eq. Ij22(l . resulting 
from the quantum fluctuations, can be neglected com- 
pared to the first term. Since 3 "'° +1 is of order unity 
around wq = 1, our previous treatment is only valid pro- 
vided that \w {ti) - 1| > 1/N. 

For smaller initial molecular population, Eq. i|16|) 
should be replaced by 



w 
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N 



3w 



w 



4{w + I) 3 / 2 



Hence, in the vicinity of Wq = 
in the w direction vanishes as 



N(w + l)V 2 (w - l) 2 
(23) 

1 the eigenvector velocity 



w = (VNa/gVS) (w - if 



(24) 



In contrast to the nonvanishing eigenvalue velocity in 
Eq. (HHJ). Substituting w from Eq. into Eq. (|T7|> . 
we obtain the action-angle, 



32 5 ^ 



n 



K-l)- 1/2 . 



(25) 



The characteristic frequency </> is now proportional to 



31, 3 



1 instead of Eq. 



(EHJ) so that AI 



F oc a . 



and 



(26) 




2 3 
rescaled time 

FIG. 8: (Color online) Atomic population fraction versus 
time, starting with a gas of fermion atoms using exact N — 
10, 100, 1000 particle calculations (solid lines) and the mean- 
field theory of Eqs. I21H (dashed lines), with 5 = 0. 



Equations (|26fl and (|2(Jfl constitute the main results 
of this work. We predict that the remnant atomic frac- 
tion in adiabatic Feshbach sweep experiments scales as a 
power-law with sweep rate due to the curve crossing in 
the nonlinear case. When the system is allowed to go near 
the critical point (i.e., when 1 — wo(ti) ^ 1/N) quan- 
tum fluctuations are the major source of non-adiabtic 
corrections, leading to a linear dependence of the rem- 
nent atomic fraction on the sweep rate. We note that 
a similar linear dependence was predicted for adiabatic 
passage from bosonic atoms into a molecular BEC |3lj . 
When the initial state is such that it has already a large 
molecular population (i.e. for 1 — wo(ti) ^> 1/N) and 
fluctuations can be neglected, we obtain a cubic-root de- 
pendence of the the final atomic fraction on sweep rate. 



VI. NUMERICAL MANY-BODY RESULTS, 
AND COMPARISON TO EXPERIMENT 

In order to confirm the predictions of Section HVI and 
Section |Vj we carried out exeat many-body numerical 
calculation for particle numbers in the range 2 < N < 
800, by Fock-space representation of the operators Ji 
and direct propagation of the man y-b ody equations |0J, 
according to the methodology of |24|. Fig. [5] shows T 
versus dimcnsionless inverse sweep rate g 2 N/a. The ex- 
act calculations are compared with a mean-field curve 
(solid green line) , computed numerically from the mean- 
field equations J2U- The log- log plot highlights the 
mean-field power-law dependence, obtained in the slow 
ramp regime a < g 2 N, whereas the log-linear insert plot 
demonstrates exponential behavior. For a single pair of 
particles, N — 2, the quantum association problem is 
formally identical to the linear LZ paradigm, leading to 
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♦ 


N=200 


A 


N=400 


□ 


N=600 


• 


N=800 




MFT 



OX) 

o 




logJg'N/a) 

FIG. 9: (Color online) Many-body calculations for the frac- 
tion of remnant atoms, V, versus dimensionless inverse sweep 
rate for various particle numbers in the range N = 2 to 800. 
The many-body results for large number of particles converge 
to the mean-field results (solid green line) computed numeri- 
cally from the mean- field limit of Eqs. 



an exponential dependence of T on sweep rate (see in- 
sert of Fig. nj. However, as the number of particles in- 
creases, many-body effects come into play, and there is 
a smooth transition to a power-law behavior in the slow 
ramp regime a < g 2 N. We note that this is precisely 
the re gim e where Eq. (|14|l can be used to estimate AI 
and F [18j. The many-body calculations converge to the 
mean- field limit, corresponding to a linear dependence of 
T on a, as predicted in Eq. 

The results shown in Fig. [3] prove the convergence of 
many-body calculations to the mean-field theory used as 
a basis to our analysis in previous sections. Having es- 
tablished the validity of this classical field theory, and 
numerically confirmed the appearance of power-law be- 
havior, we return to the experimental results of Ref. 
shown in Fig. ^ Comparison of our mean-field numeri- 
cal calculation with the experimental data (red squares 
in Fig.QJ clearly shows good agreement. However, since 
an equally good exponential fit can be found , as shown 
in the insert of Fig. ^ (dashed line), current experimen- 
tal data does not serve to determine which of the alter- 
native theories is more appropriate. We have obtained 
similar agreement with the experimental data of Ref. Q , 
but data scatter and error bars are again too large to 
conclusively resolve power-laws from exponentials. Ad- 
ditional precise experimental data over a wider range of 
slow ramp sweep rates and different particle numbers will 
be required to verify or to refute our theory. 



VII. SUMMARY AND DISCUSSION 

We have shown that nonlinear effects can play a sig- 
nificant role in the atom-molecule conversion process for 
degenerate fermionic atomic gases. In linear LZ theory, 
the precession of the two-state Bloch vector about the 
adiabatic cigenstate never stops (all the poles of the the 
integrand in Eq. I|14l) lie off the real axis), leading to 
exponential dependence on sweep rate. The nonlinear 
nature of the reduced mean-field dynamics in the large 
N limit of the many-body system, introduces dynami- 
cal instabilities. For the fermionic association case, it is 
the all-atoms mode that becomes unstable. The period 
of the precession about this mode diverges, leading to 
real singularities of the integrand in l|14|) . Consequently, 
power-law nonadiabaticity is obtained. 

While the experimental data was originally fit with LZ 
exponential behavior Q, it fits a power-law dependence 
just as well. Future experimental work with a larger 
range of sweep rates, should serve to determine which 
fit is best at low temperatures. 

The modification of LZ behavior into a T oc a 3 / 4 
power-law dependence, has been predicted for linearly- 
coupled, interacting Bose-Josephson systems [2^, l3fjj . 
Here, and in Ref. [13, we applied the theoretical tech- 
nique of |29|| . adapting it to the case of a non-spherical 
two dimensional phase space surface. The exact power- 
law for the atom-molecule system was shown to depend 
on the role quantum noise plays in the conversion process. 
When it is negligible, we find that T cx a 1 / 3 . However, 
starting from a purely atomic gas, quantum fluctuations 
dominate the dynamics, resulting in a T cx a power law. 
The same linear dependence was previously found for the 
bosonic photoassociation problem [3lj . 

Our numerical results support our analytical predic- 
tions, demonstrating how exponential LZ behavior, ap- 
plicable to two atoms, is transformed into a power law 
dependence as the number of atoms increases (see Fig.[!|Jl. 
The analysis based on Eq. i|14fl makes the differences be- 
tween the two cases transparent, relating them to differ- 
ent types of singularities. 

We note that the same power laws of T oc a 1 / 3 and 
r cx a appear in recent theoretical studies of dynamical 
projection onto Feshbach molecules [32I l3^ | . The power- 
law in |32| results from the nature of the projected pairs, 
correlated pairs giving a T oc a 1 / 3 power law whereas 
uncorrelated pairs give linear dependence, due to their 
respective overlaps with the molecular state. For compar- 
ison, the quantum noise term in our analysis corresponds 
to initial uncorrelated spontaneous emission, leading to 
r oc a linear behavior, whereas for a larger initial molec- 
ular population, correlations between emitted pairs are 
established and emission becomes coherent, yielding the 
r oc a 1 / 3 power law. The approach taken in [33j is rather 
different, based on a variant of the Wiencr-Hopf method, 
yet it results in precisely the same power-laws. 
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